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Abstract 

In this paper, we introduce a modified test statistic, by applying moving 
average method and present a new cdf estimator to estimate the joint entropy 
of the type-II censored data. We also establish a goodness of fit test statistic 
based on the Kullback-Leibler information for the type-II censored data, and 
compare its performance with the leading test statistic. A Monte Carlo simu- 
lation study shows that the proposed test statistic has greater power than the 
leading test statistic against the alternatives with monotone increasing hazard 
functions. 
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1 Introduction 

Suppose that a random variable X has a distribution function F{x), with a contin- 
uous density function f{x). The differential entropy H{X) of the random variable 
X is defined by [22] to be 



/oo 
f{x)\nf{x)dx. 
'OO 



The nonparametric estimation of H{X) has been discussed by many authors includ- 
ing Vasicek [25], Theil [2l], Dudewicz and van der Meulen [6], Ebrahimi et al. [8] 
and Bowman [4]. While Vasicek [25] and Ebrahimi et al. [^ directly obtained the 
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nonparametric entropy estimators, Theil |24j . Dudewicz and van der Meulen [6] and 
Bowman ^ derived the nonparametric entropy estimators from the nonparametric 
distribution functions. 

The entropy difference H{X) — H{Y) has been considered in estabhshing goodness 
of fit tests for the class of the maximum entropy distributions [7], [12j . 

The Kullback-Leibler (KL) information in favor of X against Y with the density 
functions g{x) and f{x) is defined to be 



J-oo 



Because I{X;Y) has the property that I{X;Y) > 0, and the equaUty holds if and 
only if / = (7, the estimate of the KL information has been also considered as a 
goodness of fit test statistic by some authors including [1], [9]. It has been shown 
in the aforementioned papers that the test statistics based on the KL information 
perform very well for exponentiality [9], and s-normality [25] in terms of powers 
compared with some leading test statistics for complete samples. 

From the censored data point of view, some authors studied the problem of 
goodness of fit test and discussed some test statistics based on a censored data 
sample. In the case of type-II censored data, if • • • , X(^r:n) be the first r order 

statistics of a random sample of size n from a distribution with cdf F{x), then the 
joint entropy of these data is defined as: 

/oo rX2:n 
■■■ /l---r:nln/l...r:nrfa^{l:n)'^2;(,r.„), (1.1) 

-oo J ~oo 

where fi-r-n is the joint pdf of • • • , X^^.^y 

However Hi...r:n, is an r multiple integral, and we need to simplify the multiple 
integral. The simple calculation of the entropy of single and consecutive order 
statistics has been studied in [19], [26] , 

[18] showed that the multiple integral Hi...r:n can be simplified to a single integral 

as: 

Hi...r:n = -(ln?i H h ln(n - r + 1)) + nHi...r;n, 



where 



r f°° 

Hi...r:n = / {1 - Fr:n-l{x)) f {x) In h{x)dx, 

n J-oo 



where h is the corresponding hazard function. 
He also showed that: 



Hl-r:n = -E (^^ ' In (J-p-^p)^ dp^ - E{{1 - Ur:n-l) Hi - Ur:n-l)), 
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and proposed the following estimator of Hi...^.n as: 

1 ^ n r r 

Hm,n,r = 'Y.^^ " "))) " " ^ " D ' 

where m is a positive integer less than or equal to | , which is called window size and 
X[i-m:n) = 2;(i:n) for i < 771 and = for i > r — m. Also [27] suggested 

another estimator of Hi...r:n based on the new cdf estimator. 

In the case of progressively censored data, first Balakrishnan et al. [3] studied the 
testing exponentiality based on KL information with progressively type-II censored 
data. Goodness of fit test based on KL information for progressively type-II censored 
data can be found in Habibi Rad et al. |13j . 

The rest of the article is arranged as follows: In Section 2, we provide a new 
test statistic: first, we explain the simple moving average method and then, we 
present our test statistic. The statistical properties of the proposed test statistic are 
obtained in section 3. Implementation of the testing exponentiality based on type-II 
censored data is evaluated in Section 4 and in this section it is shown that, for some 
types of alternatives, the proposed test achieves higher power than the leading test 
statistic. 



2 Presentation of new test statistic 
2.1 Simple moving average method 

A simple moving average (SMA) method is the unweighted mean of the previous n 
datum points. The SMA method are common tools in technical analysis |16] and 

m- 

Definition 2.1 Given a sequence ai, i = 1 ■ ■ ■ N , an n-moving average is a new 
sequence Si, i = l- -- N — n + 1 defined from the by taking the arithmetic mean 
of subsequences of n terms: 

^ i+n—l 
Si = — y Uj, 

So the sequences Sn giving n-moving averages are 

52 = ^(oi + 02,02 + as, • • • ,an-i + an), 

53 = ^{ai + a2 + as, 02 + as + a4, • • • , a„_2 + a„_i + a^), 

and so on. 
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The order selected n depends on the type of movement of interest, such as short, 
intermediate, or long term. An SMA can also be disproportionately influenced by 
old datum points dropping out or new data coming in. 

One characteristic of the SMA is that if the data have an uneven path, then 
applying the SMA will eliminate abrupt variation and cause the smooth path. In 
the next subsection, we use the definition and this characteristics of SMA method 
and present a new test statistic. 

2.2 New test statistic 

Some non-parametric estimators of (11.11) have been proposed by [18] and |27| . |27j 
expressed (jl.ip in the form Hi...r:n = — In ^^"'^^| + nHi...r:n, where 

Hl...r:n = ( ln( )dp] - E {{l - U(^r:n~l)) " f/(r:n-l))) , 

with approximating Hi...r:n by 

ln(^^_M)dp - (1 - -) ln(l - -), (2.2) 
dp n n 

in which ^ = E{U^j..n-i)), they provided its estimator as 

^ -^n(3;(i+m:n) ) ~ -^n(2^(i— m:n) ) 




H* = - > : In I ^Z^^^"^^-; . - (1 - -) HI - -) (2.3) 



n n 



where 

r{n+l)\ r-l 3;(i+l:n) - 2;(i_l:n) / 

in which m is a positive integer less than or equal to |, which is called window size 
and a;(i_„,,„) = X(i.„) for i < m and = for i > r - m. 

F~^[p) as a function of quantiles in (|2.2p is the sample path of order statistics, 
but it usually is not smooth, so relation (j2.3p can not provided a good estimator for 
(j2.2p . Then we propose to imply SMA method of proper order, say k, to smooth this 
sample path and re-consider (j2.3p for the new variables yi, • • • ,yr- So by considering 
such smoothing method, we introduce our new test statistic. 

For a null distribution function f^{x] 6), the KL distance for the type-II censored 
data is: 

h-r:nifj') = T 0) In ^^J^^^^^dx. 
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Then the KL distance can be approximated with 

r 

h-r:n{f, f) = -nHi-r-.n " In f{xf, 9) - [u - r) hi(l - F\xr; 9)). 

i=l 

Thus the test statistic based on Ii-.r-.nif'i f^)/^ and yi, - ■ ■ ,yr can be written as: 

Rm,n,r = -H^.,n,r + ^ + (n - r)y(^.„)^ ) ^) ' 

where 9 = ^ (X^I=i ?/{j:n) + (f^ — f)y[r:n)) is maximum hkehhood estimator of 9 and 

H+ = - V hi ^(»+m:n) - y(i-m:n) _ (l _ L\ ( I - L\ (2 5) 

Fn(y(i+,„:n))-^n(2/(™)) ^ ™^ ^ ' 

and 

r(n+l)V '"-I y(i+l:n) - ?/{j-l:n)/ 

in which 

_ ^(l:r)H ha:(fc:r) 

yi — fc ) 

_ ^(2:r)-\ |-3^(fc + l:r) 

2/2 — , 

-^(r — fc + l:r)H l~-^(r:r) 

Vr-k+l — fcTri , f2 Q) 

_ ^{r~k + 2:r)-\ \-^(r:r) ^ ' ' 

yr-k+2 — fcT2 > 

'^(r — l:r)"t"-^(r:r) 

yr—1 — 2 ' 

J/r = ^(r:r)- 

For y < is less than ^ and for y > yt^r-.n), Fn{y) is more than 

In the followings we assume that m in (I2.4p and (I2.5P is a positive integer less than 

or equal to | + fc, which is called window size and y(j_-m:n) = y(i:n) for i < m and 

y(i+m:n) = y{r:n) for i > r - m. 

Since I is non- negative and is zero if and only if f = f^, a.e., we reject the null 
hypothesis for large values of Rm,n,r- 

Example 2.1 For the explanation of the our proposed method, we simulate 30 sam- 
ples form the exponential distribution, consider their order statistics and censor 5 
of them from right then we plot these 25 points in Figure 1. 
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Fi gur6 1; Sample path of order statistics of 30 samples from the exponential distribution. 

Smoothing of the sample path of order statistics, and using estimator Ii2.3\) for 
new variables could provide a better estimation for \2.S\) . So we smooth sample path 
of order statistics by SMA, and provide a better estimation with greater power. 
We choose SMA of order 3 and define new variables as: 

_ ^(l:25)+a:(2:25)+^(3:25) 

yi — 3 ) 

_ a:(2:25)+a:(3:25)+3:(4:25) 
J/2 — 3 , 



_ ^(23:25) +^(24:25) +^(25:25) 

i/23 3 ) 

_ ^(24:25)+^(25:25) 
2/24 — 2 ' 

2/25 = 3^(25:25)) 

and plot the smoothed path of new variables in Figure 2. 

Figure 2 shows that the new sample path is smoother than the sample path of the 
original order statistics. If we choose SMA with order 5, it defines new variables as: 
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Figure 2: Smoothed sample path of order statistics from order 3 at SMA. 



_ ^(l:25)+a:(2:25)+^(3:25)+^(4:25)+^(5:25) 

yl — 5 : 

_ a:(2:25)+^(3:25)+^(4:25)+'g(5:25)+'g(6:25) 

2/2 — 5 : 



^(21:25) +^(22:25) +^^(23:25) +^(24:25) +^^(25:25) 

2/21 — 5 , 

_ ^(22:25) +^^(23:25) +^(24:25) +^(25:25) 
y22 4 ) 

^(23:25) +^(24:25) +^(25:25) 

2/23 — 3 , 

_ ^(24:25)+^(25:25) 
2/24 — 2 ' 

2/25 = 2^(25:25) • 

So we find an smoother sample path by plotting theses new variables in Figure 3. 

This Figure shows that even though smoothing sample path of order statistics by 
SMA of order 3 is not as smooth as using SMA with order 5, but resulting powers, 
which is demonstrated in section 4, are the same up to two digit of decimals. So 
without loss of generality, we just consider SMA of order k = 3 in (2.6\) . 



3 Statistical properties of proposed test statistic 

In this section, we consider the some statistical properties of the proposed test 
statistic in (I2.4p . The following theorem which states that the scale of the random 
variable X has no effect on the can easily 
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Figure 3: Smoothed sample path of order statistics from order 5 at SMA. 

be proved by following the line of argument in [9]. 

Theorem 3.1 Let HY...^.^ md H^,,^.^ denote entropies of the distribution of con- 
tinuous random variables Y and W , respectively, and W = kY , where k > 0. Then: 

. VariH:^^^^)=E{H+l,), 

. MSE{H+Z) = E{H+lr), 
where the superscript Y and W refer to the corresponding distribution. 
Proof: It is easy to see that: 



r{n + 1) 
for i = 1, • • • , r. So we have: 



I H + ' ' 



r — 1 w 



(i+l:n) ~ ^{i-\:n) 



fy_|_ \Y ^ \ ^ 1 ^y(i+m:n) ky(^_^-^-^ 

^ i=l -^n.(%(i+m:n)) - -^n.(A;2/(i_m:n)) 



nJ \ n/ n 



8 



Theorem 3.2 nr ~ -f^i---r:n| — m probability as n, m — t- and ^ — )• 0. 
Proof: Since [23] 

f{xi) ^ f{xi) in probabiHty, 
and [20] /(yj) f{xi) in MSE ^ f{yi) f{xi) in probabihty, consequently 

f{yi) f{xi) in probability. (3.7) 

On the other hand, we have iC(j_|_m:n)) 3;(i_m:n) belong to an interval in which f[x) is 
positive and continuous, then there exists a value x[ G {x(^iJ^rn:n)iX[i-m:n)) such that 

Pn{x {i+m:n)) ~ -^n{xii—m:n)) f, /x 

z = fi^i) 

■^(i+m:n) "^(i— m;ri) 

consequently from (j3.7p 

^niy (i+m:n)) ~ ^n{y(i—m:n)) 



y{i+m:n) y{i—m:n) 



m 



if l\ -^nix (i+ni-.n)) Pn{x {i-m-,n)) . i i -Vi 

/(xj = m probability. 

X{i-\-m:n) •^(i— m:n) 

and nr ^ n r i^ probability. 

Since [27] nr ^ -f^i---r:n in probability, the proof is complete. 



4 Implementation of the testing exponent iality based 
on type-II censored data 

Suppose that we are interested in a goodness of fit test for 

i7o:/°(x) = iexp(-f), 
i7i:/0(x)/iexp(-f), 

where is unknown. 

There are lots of test statistics for exponentiality concerning uncensored data 
[2]) [H], [H], [IS] and [T7], but only some of them can be extended to the censored 
data [5], [21], [18] and [27j. About the exponentiality test, against the one class of 
alternatives [18] showed that the power of his test statistic is greater than [5] and 
[21] also [27| showed that the power of their test statistic is greater than [18], So, we 
compare the power of our test statistic with that of the test proposed by [27] based 
on the test statistic: 



A 



■m,n,r 
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Table 1: Values of the window size m 



r 


m 


4-5 


5 


6-7 


6 


8-9 


7 


10-11 


8 


12-13 


9 


r-(r+l) 


m* (for even r) 



where „ ^. is explained in (|2.2p . 

Because the sampling distributions of Am,n,r and Rm,n,r are intractable, we de- 
termine the percentage points using 10000 Monte Carlo samples from an exponential 
distribution. In determining the window size m which depend on n, , r and the a, 
we define the optimal window size m to be one which gives minimum critical points 
in the sense of [9]. However, we find from the simulated percentage points, the 
optimal window size m. In view of these results, our recommended values of m for 
different T are listed in Table 1, where tti* — r/2 -t- 3. The critical values of Rra^n^r 
corresponding to the optimum values of m, are given in Table 2. 

Because the proposed test statistic essentially is related to the hazard function, 
we consider the alternatives according to the type of hazard functions as follows: 

• Monotone decreasing hazard: Chi-square with degree of freedom 1 (Al), 
Gamma with shape parameter 0.5 (A2), Weibull with shape parameter 0.5 
and 0.8 (A3 and A4 respectively). 

• Monotone increasing hazard: Uniform (Bl), Weibull with shape parameter 2 
(B2), Gamma with shape parameter 1.5, 2 (B3, B4 respectively), Chi-square 
with degree of freedom 3, 4 (B5, B6 respectively), Beta with shape parameters 
1 and 2, 2 and 1 (B7, B8 respectively). 

• Non-monotone hazard: Log normal with shape parameter 0.6, 1.0, 1.2 (CI, 
C2, C3 respectively). Beta with shape parameters 0.5 and 1.0 (C4). 

We consider here the sample size to be 30, and draw conclusions. We made 10000 
Monte Carlo simulations for n = 30 to estimate the powers of our proposed test 
statistic and that of the competing test statistic, for a = 0.1. The simulation 
results are summarized in Figures 4-6 and Tables 3-5. In these tables Rm,n,r and 
Am,n,r denote the test statistics of the proposed test and the test proposed by [27] 
respectively. 
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Table 2: The critical values of R, 



n 


r 


a = 0.1 


a = 0.05 


a = 0.025 


10 


5 


0.3445 


0.4253 


0.5087 




6 


0.3251 


0.4128 


0.5026 




7 


0.3136 


0.4099 


0.4929 




8 


0.3104 


0.4046 


0.4915 




9 


0.3101 


0.4038 


0.4902 


20 


10 


0.1474 


0.1913 


0.2310 




11 


0.1426 


0.1830 


0.2229 




12 


0.1408 


0.1828 


0.2197 




13 


0.1369 


0.1805 


0.2181 




14 


0.1322 


0.1773 


0.2168 




15 


0.1294 


0.1740 


0.2160 




16 


0.1281 


0.1742 


0.2161 




17 


0.1280 


0.1739 


0.2073 




18 


0.1268 


0.1649 


0.2072 




19 


0.1200 


0.1620 


0.1988 


30 


15 


0.0865 


0.1168 


0.1500 




16 


0.0859 


0.1152 


0.1430 




17 


0.0843 


0.1124 


0.1383 




18 


0.0829 


0.1090 


0.1380 




19 


0.0824 


0.1083 


0.1355 




20 


0.0815 


0.1079 


0.1311 




21 


0.0806 


0.1043 


0.1294 




22 


0.0777 


0.1038 


0.1281 




23 


0.0759 


0.1027 


0.1280 




24 


0.0741 


0.0988 


0.1275 




25 


0.0699 


0.0976 


0.1265 




26 


0.0662 


0.0977 


0.1219 




27 


0.0635 


0.0910 


0.1200 
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Table 3: Power comparison (distributions having decreasing hazard function) 



Distributions 


Al 


A2 


A3 


A4 


r 


A 


-Rui .n.v 


A 




4 

771, 71,7"' 




A 


R'TTl .Tl.T 


15 


0.0202 


0.0018 


0.0209 


0.0019 


0.0562 


0.0005 


0.0291 


0.0221 


16 


0.0092 


0.0021 


0.0089 


0.0018 


0.0313 


0.0002 


0.0246 


0.0210 


17 


0.0122 


0.0016 


0.0126 


0.0007 


0.0462 


0.0001 


0.0187 


0.0194 


18 


0.0064 


0.0009 


0.0067 


0.0013 


0.0269 


0.0001 


0.0190 


0.0184 


19 


0.0078 


0.0006 


0.0071 


0.0010 


0.0408 


0.0001 


0.0185 


0.0167 


20 


0.0023 


0.0007 


0.0029 


0.0008 


0.0189 


0.0001 


0.0169 


0.0152 


21 


0.0037 


0.0008 


0.0036 


0.0010 


0.0250 


0.0001 


0.0157 


0.0151 


22 


0.0014 


0.0008 


0.0013 


0.0006 


0.0118 


0.0001 


0.0147 


0.0127 


23 


0.0012 


0.0009 


0.0013 


0.0006 


0.0189 


0.0000 


0.0138 


0.0126 


24 


0.0006 


0.0009 


0.0008 


0.0006 


0.0047 


0.0000 


0.0092 


0.0106 


25 


0.0007 


0.0005 


0.0008 


0.0004 


0.0093 


0.0001 


0.0101 


0.0111 


26 


0.0005 


0.0005 


0.0007 


0.0005 


0.0019 


0.0000 


0.1011 


0.0093 


27 


0.0008 


0.0005 


0.0008 


0.0007 


0.0020 


0.0001 


0.1186 


0.0086 
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Table 4; Power comparison (distributions having increasing hazard function) 



Distributions 


Bl 


B2 


B3 


B4 


J" 


A 




A 




4 




A 




1 1^ 

-LiJ 


0.2190 


0.2622 


0.8890 


0.9217 


0.3974 


0.4347 


0.7065 


0.7440 


ID 


0.2383 


0.2848 


0.9118 


0.9420 


0.4106 


0.4516 


0.7254 


0.7653 


1 7 


0.2633 


0.3222 


0.9306 


0.9534 


0.4318 


0.4840 


0.7510 


0.8011 


1 S 


0.3148 


0.3504 


0.9533 


0.9689 


0.4537 


0.4958 


0.7774 


0.8278 


1 Q 


0.3369 


0.3905 


0.9645 


0.9762 


0.4721 


0.4968 


0.8034 


0.8452 


on 


0.3867 


0.4491 


0.9766 


0.9831 


0.4898 


0.5046 


0.8272 


0.8529 


91 


0.4418 


0.4915 


0.9790 


0.9889 


0.4955 


0.5201 


0.8380 


0.8581 


99 


0.5002 


0.5600 


0.9880 


0.9912 


0.5127 


0.5402 


0.8573 


0.8638 


9"^ 
Zo 


0.5799 


0.6234 


0.9927 


0.9950 


0.5394 


0.5518 


0.8690 


0.8727 


9/1 


0.6313 


0.6987 


0.9936 


0.9963 


0.5327 


0.5558 


0.8688 


0.8728 


91^ 

ZO 


0.7297 


0.7698 


0.9971 


0.9974 


0.5583 


0.5672 


0.8823 


0.8870 


9fi 
ZD 


0.7894 


0.8433 


0.9980 


0.9983 


0.5503 


0.5759 


0.8961 


0.8981 


97 
Z ( 


0.8691 


0.8978 


0.9981 


0.9990 


0.5491 


0.5831 


0.8873 


0.9039 


Distributions 


B5 


B6 


B7 


B8 


r 


A 




A 




A 




A 




15 


0.3980 


0.4362 


0.7072 


0.7455 


0.1483 


0.1754 


0.9556 


0.9750 


16 


0.4113 


0.4508 


0.7242 


0.7689 


0.1505 


0.1822 


0.9680 


0.9854 


17 


0.4290 


0.4712 


0.7529 


0.7951 


0.1671 


0.1956 


0.9816 


0.9903 


18 


0.4504 


0.4887 


0.7896 


0.8133 


0.1891 


0.2120 


0.9897 


0.9955 


19 


0.4741 


0.4909 


0.7959 


0.8300 


0.1907 


0.2311 


0.9930 


0.9982 


20 


0.4943 


0.5031 


0.8184 


0.8549 


0.2152 


0.2494 


0.9976 


0.9986 


21 


0.5045 


0.5131 


0.8375 


0.8581 


0.2408 


0.2744 


0.9985 


0.9994 


22 


0.5195 


0.5293 


0.8537 


0.8602 


0.2598 


0.2962 


0.9993 


0.9999 


23 


0.5340 


0.5423 


0.8726 


0.8788 


0.2955 


0.3233 


0.9999 


1.0000 


24 


0.5328 


0.5509 


0.8746 


0.8809 


0.3313 


0.3656 


1.0000 


1.0000 


25 


0.5626 


0.5719 


0.8795 


0.8891 


0.3752 


0.4132 


1.0000 


1.0000 


26 


0.5610 


0.5691 


0.8971 


0.8999 


0.4463 


0.5693 


1.0000 


1.0000 


27 


0.5583 


0.5701 


0.8991 


0.9001 


0.4957 


0.5221 


1.0000 


1.0000 
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Figure 5: Power comparison: monotone increasing hazard alternative at 10% when the sample size is 30. 



14 



Xable 5: Power comparison (distributions having non- monotone hazard function) 



Distributions 
r 


CI 


C2 


C3 


C4 


A 




A 




4 




A 




15 


0.9990 


0.9993 


0.6505 


0.6083 


0.3451 


0.3039 


0.0142 


0.0033 


16 


0.9992 


0.9994 


0.6518 


0.6058 


0.3334 


0.2833 


0.0056 


0.0030 


17 


0.9994 


0.9994 


0.6477 


0.6006 


0.3255 


0.2763 


0.0085 


0.0029 


18 


0.9991 


0.9993 


0.6518 


0.5867 


0.3068 


0.2481 


0.0035 


0.0033 


19 


0.9994 


0.9996 


0.6452 


0.5722 


0.2971 


0.2291 


0.0055 


0.0032 


20 


0.9998 


0.9996 


0.6267 


0.5417 


0.2659 


0.2047 


0.0029 


0.0039 


21 


0.9997 


0.9999 


0.6117 


0.5105 


0.2547 


0.1829 


0.0033 


0.0038 


22 


0.9999 


0.9999 


0.5740 


0.4803 


0.2091 


0.1556 


0.0034 


0.0049 


23 


0.9999 


0.9996 


0.5672 


0.4450 


0.1846 


0.1333 


0.0036 


0.0059 


24 


0.9999 


0.9994 


0.4894 


0.3945 


0.1443 


0.1054 


0.0056 


0.0089 


25 


0.9997 


0.9991 


0.4625 


0.3605 


0.1194 


0.0897 


0.0087 


0.0109 


26 


0.9993 


0.9986 


0.3721 


0.3164 


0.0825 


0.0716 


0.0161 


0.0183 


27 


0.9995 


0.9976 


0.3233 


0.2589 


0.0710 


0.0536 


0.0227 


0.0277 
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5 Conclusion 



As it is shown by Figures 4-6 and Tables 3-5 none of the available test statistics could 
beat others against all alternatives, but it is notable that the proposed test statistic 
shows better powers than the available test statistics against the alternatives with 
monotone increasing hazard functions, which applies to many real- life applications. 
For our test statistic against the alternatives monotone increasing hazard functions, 
Figures 4-6 show that with less remaining data, the proposed test statistic has very 
better performance than the other test statistics. 

This work has the potential to be applied in the context of censored data and 
goodness of fit tests. This paper can elaborate further researches by extending such 
modifications for other censoring schemes such as progressive censoring scheme. 
Finally, this area of research can be expanded by considering other distributions 
besides the exponential distribution such as Pareto, Log normal and Weibull distri- 
butions. 
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